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ABSTRACT 


The raoid two-dimensional slewing and vibrational 
control of the unsymmet r i ca 1 flexible 5C0LE (Spacecraft 
Control Laboratory Experiment) with mul ti -bounded controls 
has been considered. Pontryagin's Maximum Principle has 
been applied to the nonlinear equations of the system to 
derive the necessary conditions for the optimal control, 
yh* resulting two-point boundary- val ue problem is then 
solved by using the quasi 1 i neari zation technique, and the 
near-minimum time is obtained bv sequentially shortening the 
slewing time until the controls are near the bang -bang type. 
The trade-off between the minimum time and the minimum 
flexible amplitude requirements has been discussed. The 
numerical results show that the responses of the nonlinear 
system are significantly different from those of the 
linearized system for rapid slewing. The SCOLE station- 


keeping closed-loop dynamics are re-examined by employing a 
slightly different method for developing the equations of 
motion in which higher order terms in the expressions fot 
the mast modal shape functions are now included. It no 
force actuators are mounted on the beam, the modal amplitude 
responses are mere easily excited than when these actuators 
are included. System responses are dependent on both the 


f o r c e a c t u a t or i oc a t i o ns as well as t h e s t a ■- ; a nd 
weighting matrix elements. a prel iminarv study on 
effect of actuator mass on tin, closed- loop dynamic 


o n t r o 1 


the 

of large 


ii 



space systems is conducted. A numerical exarnpl 
coupled two-mass two-spring system illustrates 
changes caused in the mass and stiffness rna 1 1 i c 
closed-loop system eigenvalues. In certain cas 
for redesigning control laws previously synthes 
accounting for actuator masses, is indicated. 


e based on a 
the effect of 
e s on the 
es the need 
ized, but not 
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INTRODUCTION 


I . 

The present grant, NSG-1414, Supplement 11, continues 
the research effort initiated in May 1977 and accomplished 
in the previous grant years (May 1977 - May 1988) as 
reported in Refs. 1-15* . This research has concentrated on 

the control of the orientation and the shape of very large, 
inherently flexible proposed future spacecraft systems. 
Possible future applications of such large spacecraft- 
systems (LSS) include: large scale multi -beam antenna 

communication systems; Earth observation and resoui = 
sensing systems; orbital 1y based electronic mail 
transmission; as platforms for orbital based telescope 
systems; and as in-orbit test models designed to compare the 
performance of flexible LSS systems with that predicted 
based on computer simulations and/or scale model Earth -based 
laboratory experiments. In recent years the grant research 
has focused on the orbital model of the Spacecraft Control 
Laboratory Experiment (SCOLE) first proposed by Taylor and 
Balakrishnan^® in 1983. 

The present report is divided into five chapters. 
Chapter II is based on a paper presented at the 1989 
AA3/AIAA Astrodynamics Conference and describes rapid two- 

+ References cited Tn t h i s report are listed separately at 
the end of each chapter. 
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dimensional slewing and vibration control of the 
asymmetrical 5 COLE configuration where the beam flexibility 
is included in the model- Pontryagin's maximum principle 
has been applied to the nonlinear equations of the system to 
derive the necessary conditions for the optimal control 
where the Shuttle mast, and reflector (multiple-bounded) 
controls are considered. The resulting two-point boundary 
value problem is then solved by u^ing the quasi 1 i near i zati on 
technique, and the near minimum time is obtained by 
sequentially shortening the slewing time until the controls 
are nearly of the bang-bang type. The trade-off between the 
minimum flexible amplitude and minimum slewing time are 
di scussed . 

In the next chapter (Chapter' III) a slightly different 
method for developing the equations of motion f o i t n e tuuLE 
system during s t a t i o n k e e p i n g is presented involving a m o r e 
d i r e c t a p p roa c h in ma t r i x manipulation, and i n c 1 u d i n g higher 
o r d e i"' terms in the expressions for the mast modal s n a p e 
functions. Closed — loop responses for' t h « ^ y s t ^ ni modeled b / 
t his approac h a r e c omp a red w i t h similar r ^ s p o n s e a ^ 


presents 
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A preliminary study of the effect of actuator mass on 
the design of control laws for large space systems is the 
subject of Chapter IV. A numerical example based on a 
coupled two-mass two-spring system is selected to illustrate 
the- effects of varying the masses and stiffnesses (one at a 
time) on the closed-loop eigenvalues, and to determine what 
changes should be incorporated into the control laws 
previously designed, but net accounting for actuator masses. 

Finally, Chapter V describes the main general 
conclusions together with general recommendations. At the 
end of the grant year reported here and after submission of 
our proposal for the 1989-90 grant year 17 , the thrust of 
this research has been redirected to provide more direct 
support to the new NASA Controls/Structures Interaction 

( C 3 1 ) program. 
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II. RAPID IN- PLANE MANEUVERING OF THE FLEXIBLE ORBITING SCOLE 


The rapid two-dimensional slewing and vibra- 
tional control of the unsymmetrical flexible 
SCOLE (Spacecraft Co ntrol Laboratory Experi- 
ment) with multi -bounded controls has been 
considered. Pontryagin's Maximum Principle 
has been applied to the nonlinear equations 
of the system to derive the necessary 
conditions for the optimal control. The 
resulting two-point boundary-value problem 
is then solved by using the quasilineariza- 
tion technique, and the near-minimum time 
is obtained by sequentially shortening the 
slewing time until the controls are near 
the bang-bang type. The trade-off between 
the minimum time and the minimum flexible 
amplitude requirements has been discussed. 

The numerical results show that the 
responses of the nonlinear system are 
significantly different from those of the 
linearized system for rapid slewing. 

INTRODUCTION 

The large-angle maneuvering and vibrational control 
problem of a flexible spacecraft has been the subject of 
considerable research by many authors through different 

a PP r °aches to various structural models 1 ." 4 Among them, many 
authors placed their efforts on different control strategies 
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while using rather simplified spacecraft dynamic models. A 
few investigators have considered different and yet 

complicated structural models 3 ." 4 Among all the control 
strategies used, Pontryagin's Maximum Principle is an 
important and a basic method to such a coupled nonlinear 
dynamics and control problem. Although this method usually 
produces open- loop control strategies, it has the advantage 
of being able to handle control problems of more complicated 
structures (nonlinear dynamics and control), and it may 
prove to be useful in control- structure interaction 
problems. Unfortunately, most of the applications of this 
method to the slewing problem have been restricted to some 
simplified model, for example, a central hub with two or 
four symmetrically connected beams. Numerical problems 
appear to have limited the extension of the techniques based 

on the Maximum Principle to more complex system models 3 . 

However, by considering such extensions, we may 
encounter many interesting phenomena and produce many useful 
results. In this paper, we aim at using the Maximum 
Principle for a slightly more complicated structural model, 

namely, the 2-dimensional orbiting SCOLE 5 . The complexity of 
the present problem stems from three considerations: (1) 

more nonlinear terms than before included in the dynamical 
equations; (2) more control variables used in this system; 
and (3) the rapid slewing or near-minimum time slewing which 
may produce large flexible modal amplitudes. We hope, 
through the present analysis, to reveal, to some extent, how 
the nonlinear system is different from the linearized 
system, and how some parameters, such as the slewing time, 
and the weighting elements on the controls, affect the 
responses of the system. 

This paper consists of three parts: formulation of the 
system equations by using Lagrange's formula; derivation of 
the optimal control problem which results in the two-point 
boundary-value problem (TPBVP); and simulation of slews for 
different boundary conditions and control variables. 
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REFLECTOR 



Fig. 1 Configuration of the Planar Orbiting SCOLE 


FORMULATION OF THE STATE EQUATIONS 
System Configuration 

The Shuttle-beam-reflector system discussed in this 
paper is shown in Fig. 1. The Shuttle and the reflector are 
considered to be rigid bodies.. The beam is assumed connected 
to the Shuttle at its mass center, o . In addition, the 
reflector is attached to the beam at an offset point, a^ , 
which is x away from the mass center of the reflector, o . 

r r 

Both beam ends are considered to be fixed. 

Fig. 1 shows the structure in the pitch plane, since 
our present purpose is to analyze the planar motion of the 
system. The equations of motion in this plane are also valid 
for the motion in the roll plane, except for that case the 
inertia parameters are different. 

Three coordinate systems are used in Fig. 1 : (k g ,i Q ), 

the orbit's axes; (k ,i ), the Shuttle fixed coordinates; 

3 3 
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and (k ,i ), the reflector fixed coordinates. 8 is the 
rotation angle of the Shuttle with respect to the orbit 
coordinates. The transverse displacement of the beam from 
its undeformed position is w(z,t), where z is the coordinate 
along the k axis, and t is time. If the displacement is 
assumed to be small, then, an approxmate expression for the 
rotation angle of the cross section of the beam is, 

( z,t)=<?w(z,t)/?z. 

The free vibration of this structure can be considered 
as a free-free beam ( Bernoulli-Euler type) vibration problem 
with boundary conditions including the masses and moments of 
inertia of the Shuttle and the reflector. The partial 
differential equation for this problem can be solved by 
using the separation of variable method, in which w(z,t) is 
assumed as 

00 

w(z, t)= 2 (z)n. (t) ( 1) 

i = i g 0 

where -41. (z) is the ith mode function (shape) and q (t) is 
the associated amplitude of the ith mode. The natural 
frequencies and mode shapes for the pitch and roll motions 
are listed in Ref. 5 , and will be used in this paper. 

If the first n modes of the flexible system are used 
in the formulatiom of the dynamical equations of the system, 
the expression in Eq. ( 1 ) can be rewritten as 

w(z,t)= z t. (z)tv (t)=^ T (z)i) (t) (2) 

i ■ « * 

where . . . -41^ ] , •Q = ['Q 1 • • • ri n ] T • Then, we have, 


w = d v> / d t — i|i T t) (3) 

v)' =4 ( z , t ) = ( d\|i T /dz ) n =i4i / T 19 ( 4 ) 

d 2 v/d z 2 = ( d 2 /dz 2 ) =>v " T n ( 5 ) 

Swj=4i t (z j )Sr> ■ (6) 


8 and ri are the generalized coordinates of the system. 
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Kinetic Enerj 


The kinetic energy of the system, T, consists of three 
parts, T s , , T f , representing the kinetic energy of the 

Shuttle, the beam, and the reflector, respectively, 


T=T +T.+T (7) 

s b r 

where 

T = L I 9 1 ( 8 ) 

5 2 5 

T b =*- f L p [ (w 2 + z 2 )0 2 +v l + 2 Zw0 ] dz ( 9 ) 

J 0 

T =*-I (8 +* ) 2 +‘-m |(w 2 +L 2 ) 8 2 +2L0w +w 2 

tit ' 2 r { ' r r r 

-2x^ (9 + 43 f ) [ (i.o f +D9 ) sin4> r -9» f cos-^ ] \ ( 10 ) 


where I and I are the moments of inertia of the Shuttle 

s r 

and the reflector with respect to the attatchment points, 
respectively, m^ is the mass of the reflector, L is the 

length of the beam, w =w(L,t), and cb =w 7 (L,t). 

r r 


Potential Ener< 


The elastic potential energy of the beam is 


7-^1 r (^Ji) 2 dz 
2 Jo d z 2 


where El is the constant flexural rigidity of the cross 
section of the beam. 


Generalized Forces 

The virtual work done by the controls is 

4 

& W=u 69+Eu‘8r (12) 

1 j = 2 J J 

where u is the control torque on the Shuttle, and u and u 

i 2 3 

are the actuator force vectors on the beam, and u is the 

4 

control force vector on the center of the reflector. The 
<58 , and 5r j are the associated virtual displacements. 

From Fig. 1, we have, 

a j =u j [ cos (8 +ip . )i o - sin (8 +4) . )k Q ] 

r =( z cosfl -to sin0 )k + ( z sin0 + >.«> cos9 )i , j =2 , 3 , 4. 

J ' J J 0 J j 0 
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where u. is the magnitude of u-; z. is the location of u. 

J J J J 

along the k s axis; Vj=w(Zj,t); and $ (Zj , t) . In this 
paper, z =L/3 , z =2L/3, and z =L. After substituting these 

2 3 4 

expressions into Eq. (12), and noting the expression for Sv> 
in Eq. (6), we can get, 

4 4 

5 W=[u + S u . (z . cos$ . +w . sin* . ) ] 6 8 + Z ->ji ( z , ) u ; cos 4> j 5 7 

1 : _ . J J J J J i-* J J J 


= 0 g se+Q v 8 V (13) 

where Qg and are the generalized forces associated with 9 
and i} , respectively. 


Dynamical Equations 

After substituting the expressions (2-5) into the 
kinetic energy in Eqs. (7-10) and the potential energy in 
Eq. (11), and using the following maxtrix/vector notations. 


j pw 2 dz+m f w* =g T [ J pi|n|i T dz+m r ijj(L)'ijj(L) T ]'q=’q t M 2 

f pu ! dz+m <b J =n T [M +1 V ( L)-iJj' ( L) T g 

j ~ r r r 2 r 3 

J 0 

f L E I ( ) 1 dz=g T ( f EIW’ T dz)g=g T Kg 

Jo dz 1 Jo 

J pZwdz + I r if+rn^ =g T [ J pz^dz + l f ij)' (L)+m f Li}i(L) ]=r) T m 2 

we can obtain the Lagrangian of the system, 

Z,=Lg 2 [ I +tj t M 2 g +2m r X f (xo T cos$ f -Lsin^j^ ) ] 

• a « • • 

+9 [ig T m 2 +m f x^ ($ cos^ -w f sin<j) f -Lc|) r sinf r ) ] 

+Lg T M g -m x A w sin* ■ n T Kn (14) 

where 1=1 +1 + oz ! dz+m L* is the total moment of inertia of 

5 r J o r 

the undeformed system. The Lagrange equations, 


~ — ( — “jr ) ~~~ = Qq 


<L(^)-4k=0 

dt dr) dr) 


dt d S d 9 

of the system can be obtained in the following matrix form. 
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M(ij) 


8 


= F(9 , T) , r) ) +B(r) )u 


(15) 


with 


9 


M(D ) = 


I + i 5 T M 2 t) +2m f cos* f -Lsin* r ) 

m j +M ( 75 cos* ^ -m s sin*^ 


synunetry 

M -M sin* _ 
3 s ^ r 


(16) 


-28i] T (M z rt +m 5 co scb r -M, 15 sin* r ) +m f ** ( Lcos*^ sin^^ ) 
9 1 (M^ 15 +m a cos* r -M ( i) sin* r ) + (* y -28 )i} cos* f -Kri 


(17) 


cos* 4 — sin* 4 
H,(z 4 ) cos* 4 

(18) 

where u=[u u u u ] T is the control vector. Other 

t 2 3 4 

notations used in these equations are 

M =m X Y (L)V(L), M =M +M T , M =M -M T ; 

M 4 =m r x f y (L)y T (L) ; 

m s =m f m ( L ) + Li y ( L ) ] , =m f X f hjj ( L ) - Ly' ( L ) ] 


B= 


1 z 2 cos * 2 -w 2 sin * 2 
0 y(z t )cos * 2 


z cos* -v) sin* 

3 r 3 3 r 3 

1 ||(Z 3 )cos * 3 


We need the following linearized version of Eqs. (15) 
to compare the responses of the two systems. 


Z ! ml' 


’0 " 


'o I 0 T ' 

__ 1 __L 



= - 


m ! M 


T) 


0 ! K 

2 1 3 




— 


'l z 

2 

z L 

3 

0 1) ( z 2 ) 

^(z 3 ) ip(L) 


(19) 


For convience, by introducing the notations 


y[ = [e , 15 T ] = [y 1 t 


y ik I. k=n+l. 


y =y 

2 2 2 1 



Eqs. (15) can be rewritten in the state form 
y =Y 

2 1 2 1 

y 2 =M" 1 (tj ) [ F(n , y 2 )+B(i3 )u] 


} ( 20 ) 
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DERIVATION OF THE OPTIMAL CONTROL PROBLEM 


Obiecti ve 

The purpose of this paper is to find the optimal 
controls which rapidly drive the system from an initial 
state, y(t=0), to a final required state, y(t=t f ). Since the 
magnitudes of these controls are, in pratice, bounded, the 
optimal controls for the minimum time slewing problem are 
usually of the bang-bang type. However, this kind of control 
will generally introduce large flexible amplitudes. 

Therefore, a near-minimum- time slew is of primary interest 

to us. 

Necessary Conditions 

Instead of starting from the minimum time control 
problem, we set out to deal with the optimal control problem 
with a quadratic cost function, 

y=if tf ( y* Qj y, + y[Q 2 v t +u T Ru)dt (21) 

2j o 

where Q , Q , and R are weighting matrices, t f is the given 
slewing' time. This kind of problem has been considered by a 
list of authors. However, in their analysis, t f is fixed and 
there is no limitation on the magnitude of the controls. On 
the contrary, in the present problem, the slewing time t f is 
no longer fixed, because we want to find a rapid slew or a 
near-minimum- time slew. The magnitudes of the controls, u, 
are also bounded, 

lu. |<u. , i— 1 ,2,3,4. ( 22 ) 

l l 1 i/ b 

Our strategies to solve this problem are described in 
the following. First, the necessary conditions based on Eqs . 
(20-21) are derived. Then, the costraints, Eq. (22), are 
imposed on these necessary conditions to modify the 
controls. Finally, in the solution process of the resulting 
TPBVP , the slewing time is shortened sequentially, in order 
to find the near-minimum-time slewing. As we have discussed 
in Ref. 6, when the slewing time is shortened, the optimal 
control, will approach the optimal control of the minimum 
time slewing problem, that is, becoming the bang-bang type. 
It is clear that, when the controls approach the bang-bang 
type, the value of the index J in Eq. (21) will increase and 
approach its maximum value. 
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( 23 ) 


The Hamiltonian of the system is, 

ff=^(y[Q, y t + yjQ 2 y t + u T Ru)+x[y i +\]'m' 1 (f+bu) 

where x and x are the costate vectors associated with y ( 
and y 2 , respectively. By using the Maximum Principle, the 
necessary conditions for the unresticted optimal control 
problem are the dynamical equations (20) plus the following 
differential equations for the costates. 


X =-M =-Q y -iI T M‘'\ > ( F+Bu ) - <X T M" 1 — >u 

1 <?y, 1 1 <?y, 2 *<?y ( 1 dy i 


X =-— =-Q y -X M' 1 x 

2 w 2 2 1 ir 




<?y, 


<?y. 


•jtt '~[ r 7 r :> ■ rtf 


where <xj44?-' >=fiM-X, , 4 


r 


represents a 


special matrix (similar for <\J M~ 1 >); as well as the 

expressions for the optimal control, 


(24) 

(25) 


|tf=0, u=-R' 1 B t M* 1 X 

The control rules in Eq. 
following expressions 7 , 



r -u, 

b ’ 

if 

u. <-u, 

l> C 1, 

b 

u, =- 

If 

u. 

c » 

if 

3 

V 

u 

zi* 

b ’ 


u, 

- If 

b * 

if 

_c 

n 

V 

b 


(26) 

(26) are then modified by the 
where 

. u. =- (R“ 1 B t M' 1 x ). , (27) 

C 2 i 

i=l ,2,3,4. 


By substituting the control expressions into the dynamical 
equations (20) and the costate equations (24), we can 
obtain a set of 4(n+l) differential equations for the states 
and the costates. To obtain the control, u, we need to solve 
this set of differential equations with the 4(n+l) given 
boundary conditions: y(t=0) and y(t=t f ). This problem is 
called TPBVP because the B.C.s are specified at the two 
ends of the slewing period. 


Solution of the TPBVP 


The quasilinearization algorithm and the method of 
particular solutions are used to solve this nonlinear TPBVP 9 . 
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NUMERICAL RESULTS 


Some common parameters of the SCOLE used in this paper 

are , 

EI=4xl0 7 lb- ft 2 , p =0 . 09554 slug/ft, L=130 ft, 
m =6366.46 slug, m =12.42 slug, 

s r 

u b =10,000 ft- lb, u b =u 3b =10 lb, u 4b =800 lb. 

Other different structural parameters are listed in Table 1. 

Table i 

STRUCTURAL PARAMETERS OF THE 3-D SCOLE 


Roll-Axis 

Pitch-Axis 


905,443 

6,789,100 

slug-f t : 

18,000 

9,336 

slug-ft : 

32.5 

18.75 

ft 

0.319954 

0.295016 

hz 

1.287843 

1.645292 

hz 

4.800117 

4.974182 

hz 


All the numerical tests done in this paper are 
rest- to- rest slews, that is. they use the same boundary 
coditions for the states: i)(t=0)=0, i)(t=t f )=0; 0(t=O)=O, and 

8(t=t f )=9*, where 9* is the required slewing angle, ranging 
from 20 deg to 180 deg. All these slewings can be divided 
into the following 3 groups. 

Group 1 

In this group, only the Shuttle control torque has been 
used, i.e., u=u ( . The weighting matrices Q, =Q 2 =0 and the 

weighting on u , r =10” 5 . Figs. 2 show the near-minimum-time 
slewing about the roll axis, through 90 deg (Fig. 2A) . The 
near-minimum- time , t f , has been calculated to be 27.8 sec. 
The control torque is near the bang-bang type (Fig. 2F). The 
maximum amplitude of the first mode of the linearized system 
is 9.2 ft (Fig. 2B), which is less than 10% of the total 
length of the beam. The first modal amplitude response of 
the nonlinear system has a shape similar to that for the 
linearized system, but with a shifting of the amplitude. The 
second mode and the third mode of the nonlinear system have 
quite different time histories from their linearized 


counterparts (Figs. 2C-D) . The rotation angle, $ , and the 
displacement, u> , at the reflector end of the beam are also 
plotted in Figs. 2A and 2E. They have shapes similar to the 
amplitude of mode 1, because the first mode dominates the 
deformation of the beam for this slew. 

The slewing about the pitch axis has responses similar 
to those about the roll axis. To make a comparison, the 
results of many other slewings in this group are listed in 
Table 2. r\ is the maximum value of the first modal 

'1 IHJ 

amplitude of the linearized system. Note that the number of 
vibrational cycles of the first mode increases as the slew 
, * 

angle, 8 , increases. 


Table 3 

RESULTS OF GROUP 1 



Roll 

Axis 

Pitch 

Axis 

0 *(deg) 

t f (s) 

ij (ft) 

' l ma x ' ' 

t f (s) 

i), (ft) 

* 1 m a x 

20 

15.99 

7 . 5 a 

31.85 

2 . 8 a 

45 

20.56 

9 . 5 a 

48.29 

2.8 b 

90 

27.80 

9 . 2 a 

67.05 

2 .Q r - 

180 

40.14 

9.5 b 

95.23 

2.8 d 


Q One cycle. b Two cycles with one big peak and one small 

peak. c Two cycles with two equal peaks. d Three cycles with 
two big equal peaks and one small peak. 

Group 3 

In this group, the force on the reflector, u^ , is added 
to the system. The weighting on the states, Q ( and Q 2 , are 

still chosen to be zero, and r ( =10 -s . The effect on the 
slewing responses of adding the control force may be 
analyzed by changing the values of t f and/or r^ , the 
weighting on u^ . Since the first modal amplitude dominates 
the deflection of the beam, our main concern will 
concentrate on the variation of the first modal amplitude. 

To illustrate the effect of the parameters, t f and r^ 
on the time response of the first modal amplitude, let s 
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consider a special case without lose of generality, i.e., 
the 90 deg slewing about the roll axis, the same case 
plotted in Figs. 2 but with the control u 4 added. In Fig. 

2B, the time response of the first modal amplitude can be 
approximately expressed as -tj ma x sin(2nt/t f ). This response 

is 180*out-of phase with the control u t (Fig. 2F) , because 
of the inertia effect of the flexible beam. However, when u 4 
is added to the system, the torque produced by u 4 will 
accelerate the slew and balance the deflection of the beam 
produced by u ( . 

It is not hard to imagine, from the physical point of 
view, that when u 4 increases to a large value, the response 
of the first modal amplitude may be in-phase with U ( (or 
u ), i.e., (t)sri sin(2nt/t. ). Therefore, between the 
small values and the large values of u , there must exist a 
critical value at which the phase of the first modal 
response changes from out-of-phase to in-phase. It is also 
expected that during the "phase-change" period, the maximum 
value of the first amplitude becomes a minimum. This 
conjecture, fortunatly, has been proved to be true in our 
calculations . 

One way to change the value of u 4 is to change the 
value of r 4 , for fixed slewing time t f . Another is to change 
t while mantaining r 4 fixed. These results are plotted in 
Figs. 3A-3B . We should point out that for large values of r 4 
(Fig. 3A) or large values of t f (Fig. 3B), u 4 is small and 
the response of the first modal amplitude is out-of -phase . 

On the contrary, small r 4 or t f results in large u 4 and, 
therefore, in-phase response. In each of these cases, a 
minimum value of n exists. It is also interesting to 

know that, at these critical values of r 4 or t f , -g(t) 
experiences two oscillation cycles with two equal peak 
values (or valley values) of the linearized system, i.e., 
n(t)sn sin( 4n t/t„ ) . The dotted lines in Figs. 3A-B 
represent the nonlinear system responses. The nonlinear 
response has a shift from the linear response, especially 
when t or r is reduced. Also, we have observed that, at 

f 4 

the critical points, although the amplitudes are small in 
value, the linear and the nonlinear systems have quite 
different time response histories. 
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A complete relationship between the three parameters, 
tj i m a x » t f , and r 4 can be investigated through the 
3-dimensional surface in Fig. 4. The lower ditch on this 
surface represents the minimum value area of ij . Although 
the global minimum value of rj occurs when t A is quite 

large, there exists a local minimum value, t> =0.41 ft, 

'I nan 

around the middle of the ditch, where t =23.881 sec and 

f 

r 4 =0.86. This important point can be chosen as the trade-off 
point between rapid slew and small amplitude requirements, 
at which neither t f nor is too large. The response 

shapes of the first modal amplitude for the different values 
of t f and r^ are different. In the hilltop areas, only one 
vibrational cycle of ij (t) exists, but along the deep 
valleys of the ditch, tj (t) has two vibrational cycles with 
two equal peaks. More surprisingly, at the local minimum 
point mentioned above, tj (t) experiences three vibrational 
cycles with three equal peaks. The responses for this case 
are shown in Figs. 5, where the linear and nonlinear systems 
are quite different in spite of the small fnodal amplitudes. 

Group 3 

Based on the example shown in Figs. 5, the controls, U 2 
and u are added to the control system in this group. The 
associated weightings on these controls are r 2 = 10.0 and 
r 3 =20.0. Also, the weightings on and 13 are selected as 
200.0 and 1000 . 0 , respectively, to show the further 
reduction of the modal amplitudes. These results are plotted 
in Figs. 6 . Compared with the results in Figs. 5, the modal 
amplitudes have been slightly reduced and the maximum value 
of u has been reduced due to the addition of u and u . 

4 2 3 

Note that u ; is not shown in Figs. 6 because of its 
similarity to that in Figs. 5. 


C0NCLUSI ON 

The Maximum Principle has been applied to the rapid 
slewing problem of the planar flexible orbiting SCOLE . The 
dynamical equations used contain more nonlinear terms than 
those used by other authors, and the responses indicate the 
large differences between the nonlinear and the linearized 
systems, not only in the rapid slews where large modal 
amplitudes are involved, but also in the small -amplitude 
slews. The analysis between the relationship of the 
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parameters, tj , t f , and r^ , indicates that the conflict 
between the rapid slew and the small flexible amplitude 
requirements may be compromised for multi-input control 
systems. The effects of these parameters on the 
3-dimensional SCOLE model slewing responses need to be 
investigated . 
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00 Deg Boll Axle Slew, Using Shuttle Torque Only. 
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Figs. 2 90 Deg Roll Asia Slew, Using Shuttle Torque Only, 
(continued) 
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Figs. 3 Variation of Mode 1 vs. r4, and Tf. 
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Fig. 4 Relation between First Modal Amplitude, 
Reflector Control Weight (r4), and Slewing Time. 



x iooo ( rr-LB ) 



90 Deg Ron Axle Slew, Tf-23.881 a, r4«0.86. 
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90 Deg Roll Axle Slew, Using All Controls. 







III. CONTROL OF THE ORBITING SCOLE WITH THE FIRST FOUR MODES 


A. Formulation 


In order to complete the calculation of the elements of the state and control influence matrices 
for the orbiting SCOLE system linearized about the nominal station keeping motion, we list all 
equations of the system which are based on the formulation of Ref ( 1 1 as follows: 

1. Generic Modal Equations of the Beam 


A n + < A n - 


2 <*> 


1 G 1 ( Pn ) *1 2 + “j“ G 2( Pn ) h I + f 0) () G 3^ Pn )t 1 2 


L 


r " G 2(Pn^3 + L ° 2( ^ >1f| 1 — Gl ( F 


where 


(n=l ,2,3,4) is a time dependent amplitude of the nth mode. 

^ i (i=L2,3) are angular displacements about roll, pitch and yaw axes. 

G,( P„ ) = f 3 ( P„ )A 4- f 4 ( P n )B , „ + f s ( p„ )C „, + r ( f p„ )D„, 

°2< Pn ) = f 3< P„ > A 2„ + f 4< P„ )® 2n + r 5 ( P„ > C 2„ + ^ P„ » D 2n 
G 3<Pn)= f 3<Pi>3„ + f 4<P>3n 


sin( p L) 

f 3 <Pn)= — f- 
P" 

cos ( P n L) 

P 


f 4 C P n ) = 


— + 


Lcos ( p n L) 
Pn 


Lsin ( P n L) 


Pn 


P 


Lcosh ( P L) 
f 5(Pn)=— n 


f*( Pn ) = 


Pn 

Lsinh ( P L) 


sinh ( P n L) 

P 2 

r n 


cosh ( P n L) 


Pn Pn Pn 

F = F x ( V 3x S xn(- L ) + ^x S xn (" 2L/ P) +V tx S xn I 

+ F yf V 3y|n(- L ) +V 2ySyn<- 2L/3 > +V 1y V ,73)| 


( 1 ) 
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S xn ( Z) = A ln sin( P n Z) + B ln cos ( P n Z) + C 1n sinh( (3, 7 ) + D ]n cosh ( P n Z) 

S yn ( Z) = A 2n sin( P n Z) + B 2n cos ( p n Z) + C 2n sinh( p,Z> + D 2n cosh ( p n Z) 

e n ( z ) = A 3n sin ( Pn Z) + B 3n C0S < ^n Z ) 

P = P 2 J E 

Pn Pn V G A 

2. System Equations without Flexibility and External Forces 

1 ] XX - ^ 2 r xy - ^ 3 - W 0 T I 3 ( ] xx ~ [ yy + ] 77 ) ~ “d 1 ' ? ’yz 

- i( \ 77 - l yy ) -cog T ! 3 f X7 - ti 2 I xy 

B 2 T yy +T ^l I xy + ^ 3 l^z + “o 1 ! 3 f xy "W^l’xy 

+ W 0 TT| 3 I y7 + 3 W Q T| 2 ( I xx \ y7 ) — 

^ 3^7 ~^l l xz "^2 T yz + “o 1 ! T xx _1 yy + ‘zz ] 2 ! xy 

- 4w 0 t Ii I xz +3w 0 T l2 I y7 -"> 0 ^ “I vv ) 

where 

C* = Cl +I R1 + ^ +M^L 2 +Y 2 ) 

= ',2 +, R2 +! ^ + MR( L2+X2 > 

= ' s 3 + I K3 +M R< X2+Y2) 

Cy = M R XY 

T X7 = ! s4 + M R* L 

!yz = M R YL 


= 0 ( 2 ) 


(3) 


= 0 (4) 
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3. System Equations with the First Four Flexible Modes 


^ | J xx - ^ 2*xy - ^ 3 [ xz ~ 3< ! xx “ T yy + l ?y) ~ M J v/ 

-^oV^Z* I yy)- W 0^3 I xz " 3w O T l2 , xv 

+ ,lt A n d ln -J 1 A n d 2n + j., A n d 3n 


= T. 


2 r y y +T i l ! xy + ^3 T yz " ^“oV- +T| 3 co O I xy _3( °0 V' 


+ “o t i 3 I yz + 3 (o§t, 2 (i xx -i Z7 ) + n J ( A n d 


4n 


^ , A n d 5n = 

n~ 1 


^3 1 ZZ l T xz “ ^2*yz + f xx “ ! yy + I 7z> -(0 o^2 l 

- 4t °oVxz + 3w 0 7 l2 r yz " w 0^3( ! xx " 1 


VV 

+ |, A . d 6„ + j/n"7„ 'J, A , d «n 


T. 


where 


, T , 

I , 1 

. 1 are same 

as 

in 2 

yy ’ zz 

xy ’ yz 

xz 




M x U x + 

FyL[ V , y / 

3 + 2V 2 y / 

3 + 

V 

3y 1 

M y u y + 

FxMV, x / 

3 + 2V 2x / 

3 + 

V 

3x 1 

M z U z + 

XF y V 3y 

- YF x V 3x 





B. System State Equations 


In this section we recast all system equations (1-7) into matrix form. 
Let 


T 

X = I t) h T) A A A A r) r| T) A A A \ | 

1 1 2 3 1 2 3 4 1 2 3 I 2 I 

as a state vector and 



U=(Mx My Mx M, V 3x V 3y U x U y U x I 

as a control input. We then set up the system state equations hv two different methods and get the 
state matrix and influence matrix, respectively. 


1. Method of Ref [ 1 ] 

Generic modal equations of the beam: 


V 


V 


*il 




^1 

A 2 


> 

to 




0 2 

. . 

+ [dJ 


+ [dJ 

^2 

+ [dJ 

0 2 

+ [dJ 

A 3 


A 3 

L 2 J 


L 3 J 1 



A 4_ 


A 4_ 


1 

1 


1 


0 3 


= [f] 


Vi x 
VI y 
V 2 x 
V2y 
V3x 
V3y 


System equations without flexibility and external forces: 



^1 




Ol 

[eJ 

^2 

+N 

^2 

-[eJ 

^2 


^3 


^3 


^3 


= 0 


System equations with the first four flexible modes: 





X 


T 'l 


A 1 


Tli 


+ [eJ 

A 


+ [e 4 ] 

A 2 

-[eJ 


N 

o 2 

A 3 

+ [ E J 


A 3 

02 


^3 


_ A 4_ 


0 3 


_ A l 


0 3 




1“ 

1 

r 

“1 


r "i 



•[=J 


" A i" 


VI X 


l! 

A 


Vly 


X 

A 2 

A- 

- [mJ 

V 2 x 

V 2 y 

+W 

U y 

3 

V3x 



A 4_ 


V3y 


u* 


We then recast eq(9) by inverting the matrix fell . 


( 8 ) 


(9) 


(10) 
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( 11 ) 


*1 

r “ 1 i 


r “1 i 


^2 

^3 

- -[e, e 3 ] 

^2 

^3 

- [ E i eJ 

i 

-zf J 


After substituting eq(ll) into eq(8), the result is, 


V 

^2 

= -N 

V 

a 2 

A. 

-[d 3 - D 2 E 1 E 3 ] 

1 ” 

ct 

. p- • r 

1 

-[D 4 -n,F, Ej] 

^2 

+ [f] 

Vlx 

Vly 

V2x 

V2y 









V3x 



A 4_ 


^3 


^3 


V3y 


or briefly 

A = [C J A +[C J b +[C J r, +[f] V (13) 

where 


A = 

[ 

A 1 

A 2 

A 3 

aJ t 

A = 

[ 

A 1 

A 2 

A 3 

a 4 ] T 

T1 = 

[ 

b 

1 

b 

2 

b 

3 

r 

b = 

[ 

b 

1 

b 

2 

b 

3 

f 

V = 

[ 

V 1x 

V iy 

V 2x V 2y V 3x V 

[cj 



[dJ 



[cj 



[«V 

- D 

2 E 1 E J 

[cj 

= 


L D 4- 

- D 

2 E 7' E sl 


Then eq(13) without the external Forces is substituted into eqMO) with the result 

'r\ = [cj A +[cj t, +[c f J t, + [C J A + [Mjv ' F'lju (14) 

where 
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* - M, \ r-,] T 

A = [ A, A 2 A 3 A 4 ] 

y X 

U= [U x Uy U 2 ] 

[c J [e 1 ( E 3 +E 2 C 2 )] 

[cj- - [e;'e 4 ] 

[c J = - [e 7 '( E 5 + E 2 C 3 )] 
[cj- - [e T'< E 6 +E 2 C,)] 


[mJ = [e;'m,] 

[mJ = [e;'m 2 ] 

Eqs(14) and (13) may be combined as follows: 


T) 


1 

n 

4^ 

c 5 ' 


b 

+ 

_C 6 

C 7 


h 

+ 

M, 

m 4 


V 

_ A_ 


1 

n 

0 _ 


B A_ 


. C 3 

C l. 


A _ 


F 

0 _ 


u 


The system state equation becomes 
X = [a ] X + [ B ] u 

where 


0 

0 

I 

0 


n 

0 

0 

0 

0 

I 

and [ B ] = 

0 

0 

C 6 

c_ 

C 4 

C 5 

VI 

M 4 

C 3 

c . 

C 2 

0 


F 

0 


\ 


(15) 


2. Direct Method 

The generic modal equation (eq( 8 )) and system equation (oqMO)) may be directly combined 
to yield: 
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Eq(16) may be rewritten, following the inversion of the acceleration coefficient matrix, as 



where 



We can get the system state equation from eq(18), that is 
X = [ A ] X + [ B ] U 

where 



( 19 ) 
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C . Control Synthesi s - LQR 

The svstem state equation can be represented as 

X = (A] X + C B ] U < 20 > 

An LQR cost function is selected as follows: 

•J = / (X T QX U t RU) dt 


( 21 ) 


The optimal control, U, based on the LwR theory is 
given by 

u=-[R''b t p]x < 2 2 j 

where P is the positive definite solution of the o t s a d v 
state Ricatti matrix equation: 

P A + A t P * P8R' 1 '3^P f Q 3 0 v33) 

The closed loop system equation becomes 

X 3 [A-BK] X (24) 

Let X ( 0 ) be an initial state vector. Based on some 
assumed Q and R penalty matrices, the closed loop clvnami 
responses csin b 6 s i ni u i 3 *t 6 ci 3 3 
X ( t ) = et A_eK l t XiO) 

which i s based on the feedback control 9 i v « n o v 
ij ( t ) ~ - !< X ( t ') 

The total tornue impulse about the three axes an 


< 25 ) 


(26) 


oo 

T < ’ 

IMt) 1 

dt 

'V z /" 

|Ty(t) | 


o 



T z = /°° 

| T z ( t ) | 

1 dt 


( 29 ) 


w h 3 l 


'T;<i tf 
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— — 
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1 C 
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Numerical Results 


D . 

The OR AC L 5 control software in the IBM computer' system 
was used to calculate the state matrix [A] and influence 
matrix [B] and to simulate the closed loop system responses 
as well as the total torque of the system for given sets of 
initial conditions. 

We select the force factors, F x = F y = 1 and torque 

factors M x = M y = M z = 1, which means that the components 
of V and U reflect the actual actuator force and Shuttle 
torque values. According to the SCOLE configuration and 
parameter values (listed in the Appendix), the [A] and IB] 
matrix values of the method of Ref. [1] and the direct 
method are listed in Tables 1,2,3, and 4. 

We select the initial states 

X | r 0 ) = 6 degrees 

x 2 ( 0) = X 3 V 0 ) = X 4 (0; = X 5 1. 0 ) * X 6 ( 0) - X 7 ' 0 ) = 0 

and the diagonal weighting matrices as: 

trace Q - [ I 0 ^ > 1 0 ? , I 0 ^ , 5 x I 0 ^ . 1 0 ^ . 5 ,< 1 0 ° . 5 x I 0 - • 

10,10,10, 10,10,10.10] 

trace R - (100,100,100,100,100,100,0.001,0.001,0.001] 

The simulation of t h - optimal closed- loop system 
responses, using both the method of Ref. I. i j and the direct 
method . are p lotted in rigs. 1 . 2 , 3 and k . 

The tot a i control torque - 1 mo u ! se s c r the system 
rs. :: u i : ■ e d abut:- the roll -axis air Ou.528 bv the 

method of Ref. I 1 ] and 3 : ,2 50 ft -lb -sec r-.r th 3 direct 

r.-:. : 

OF POOR ; o .. o 
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method. The torques needed about the other two axes are 
much less than the components about the roll axis. Also, 
the maximum torques of the system are 6,525 ft-lb for the 
method of Ref. [I] and 5,302 ft-lb in the direct method. 


E . Concl usi ons 

1. Bv comparing the results of the method of Ref. [1] and 
the direct method, it is seen that the results ars similar 
to each other. 

2. In the responses resulting from the direct method for 
the same initial displacement about the roll axis, it is 
seen that the first four flexible modes are general ly 
excited more than for the results of the method of Ref. [1]. 

3 . If no force actuators are added to the beam and 
reflector complete damping of the modal responses requires a 
much longer time (Fig. 6) than when the force .actuators are 
utilized together with the Shuttle torquers (Figs. 2 and 4). 
However', the use of force actuators results in initially 
larger overshoots as compared with the case depicted in Fig. 
6 . 

U . The system responses are dependent o n the force actuate r 
I o cation s a nd 1 1 *. a w e i g h t i n g m a t r i ces ( Q . R ) v a lues. Suitable 
values of the penalty matrices and actuator locations 
should be selected so that the system control becomes 
opt l ma : . 

3. From the system .analysis, we rind the f‘ sxibilitv of the 
3 C 0 L E s v s t &m i s n c t g , : i t I v a v: c i t e d d u r i r. y t v p i c a l s t a t i o n 
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System responses and the total torque 


keeping operations, 
impulses needed are similar to the rigidized SCOLE system 
( see Ref . [ 1 ] ) . 
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NASA Grant NSG-1M4, Suppl. 9, January 1938. 
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Transient Responses about Roll, Pitch, and Yaw axes 
(Method of Ref. (1) ) 
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Fig. 2 : Transient responses of the First Four 

Flexible Modes (Method of Ref.(l)). 
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Fig. 3 : Transient Responses about Roll, Pitch and Yaw Axes 

(Direct Method) 



Transient Responses of the First Four 
Flexible Modes (Direct Method) 






Pig. 5 : Transient Responses about Roll, Pitch and Yaw Axes 
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Flexible Modes without Force Actuators 




DRAWING OF THE SCOLE 




CONFIGURATION 
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APPENDIX 


Format of Submatrices 




- 3 “ > 0 I xy 

-3»Slxy 



-(On I 
0 X 2 


“> 0^2 


L ■ 4 "o i „ 



3»§Iy, 

r d =' 

d 3 2 

d 33 

d 3-4 

° 

0 

0 

0 

L d «. 

d 82 

d S3 

d 8 < _ 


( XX ^uv) 


yy - 
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0 

0 

0 


0 ” 

0 

0 


CO' 



g 2 ( 0 ,)/l -G,(e,)/L 

G 2 (o 2 )/L -g,( 0 2 )/l 

g 2 ( 0 3 )/l - 0 ,( 83 )^ 

G 2 ( 0 4 )/L “Gj ( 0 ^ )/L 


0 

0 

0 

0 



0 

0 

0 

0 


2co 0 G 3 (o' )/L 
2co 0 G 3 (0 2 )/L 
2<o 0 G 3 (o 3 )/L 
2co q G 3 ( 0 ^ )/L 


-2co q G 2 (o, )/L 
-2co 0 G 2 ( 0 2 )/L 
-2co 0 G 2 (o 3 )/L 
-2co 0 G 2 (o 4 )/L 



4<o2G 2 (0. )/L 
4co2G 2 (o 2 )/L 
4co 2G 2 (0 3 )/L 
4co2G 2 ( 0al )/L 


— 3co^G j ( 0 j )/L 
-3<ogG,(0 2 )/L 
-3co2G 1 (o 3 )/L 
-3^,(0^)^ 


0 n 


0 

0 

0 


- F , S<i (-L/3) 

F x S x 2 (- l /3) 

F . s ,,<- L/3 > 

F , S x,<- L / 3 > 


F y S y ^ ( ~L /3 ) 

F A 

F yS y V L /3) 

F x S x. 

FySY^-V 3 ) 

F x S x 

F y S y ( -L/3 ) 

F x S x 


F ySy ( -2L/3 ) F.S^f-L) 

Fy S '(-2 L/3 > F,S, 2 (-L) 

FyS^I-aL/ 3 ) F.S.yl-L) 

FyF y 3 ( -2L / 3 ) F.S^I-L) 


(-2L/3) 

( -2L/3 ) 

( -2L/3 ) 
(-2L/3) 

I 

Vy,<- L) “ 

F y\<- L > 
F y S (-L) 

F y s y 4 *“ L) J 
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0 


0 


0 



F y L/3 
-F x L/3 0 

0 0 


F y 2L/3 


-F x 2L/3 0 

0 0 




0 


0 


0 

0 




2. System Flexible Mode Shapes 
(1) Method of Ref.[l] Equation (For nth mode) 
d ln =M R [ “ LS ny ( - L ) -XY® n ( - L ) ] +Mf 2 ( B n ) /L 

d2n= M »“ol YS n.(- t ‘> 1 
d 3n = w§(Mf 2 (B n )/L-M R Xe n (-L) ] 

d 1n =M£, ( ) /L+ [ I R; ‘Mj ( x2 + L 2 ) ] ( - L ) -MjLS nx ( -L ) 

d : n ="o M R X|L ®B<- L > + 2 S nx<- L >> 

d S n =M R[ XS n y(- L )- YS n » ( - L)1 

d 7 =«V> II XI. B (-t.) 

d sn="^S n ,(-L) 



(2) Direct Method's Equation (For nth mode) 


d ,n =Mf 2 (»„ > / L *V S ny < ' L > *«R XL *n < ‘ L > + V^ ' „ y < L > 
- M R XyS '„,(- L >* I « 1 S 'yn(- L > 

d 2n = ^ w oH R YS nx (-L) + (-2w 0 M% R Y 2 -to 0 I Ri +co 0 I R3 -w 0 M R X 2 -«o 0 I R3 )* n (-L) 

d 3n=“5 1 Hf 3 < “n ) / L tM R XYS ' < - L > - V JS ' < - L > 

♦*^<- L > + VVVV- L)l 

d 3n= Mf , ( °n > / L -"R LS n, ( - L ) * J R 2 S ' ( - L > + V L ®n ( ^ > 

-«R XYS 'ny(- [ ‘) t V 2S '„,<- L > 

d Sn = - M R"0 X S„ x ( - L ) +M,« 0 YLS ' „ y ( - L ) -M^XLS ' n> ( - L ) +M R « 0 XY®„ ( - L ) 

<WV®Ry < - L > - V s ** < -I !• ) + V Y ®» < - L ) + V 2 »n < ‘ L > + V" ( ' L ’ 

d 7„=' Vo Y2s ' < - L ) L !r 2 S ( - L > + “o Ir 3 S • „ y ( - L ) 

M R XLa, 0 e n (-L)+<o 0 I R] S' ny (-L) 

d 8n = M R"i YS ™<- L > + VS XS ny(- L )-“? I R | ® n <- L >*“i I R 2 9 n(- L ) 
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where 


■■<°n> =A .n( 


Lcosb„L sine L 

H n 


B n B 


sinha L LcoshB„L 
+C 1 —~ 


l0 n L 1 + R r L3in B n L , COSB n L ^ 1 >| 

2 J ln [ °n ®n y 

ihB n L 1 +n r LsinhB n L coshB n L ^ 1 ^ 

n J ln [ B n 02 b 2 J 


r Lcosb„L sins„L^ „ f Lsino L coso L 1 ^ 


+c 


2n 


( 


sinhfl n L Lcoshs„L 


ihB_L 'N f 


LsinhB„L coshB n L i 
n n 


0 


S 'nx (_L)=B n [A ln COSB n L+B ln sinB n L+C ln COShB n L+D ln sinhB n Ll 

S ' ny ( -L ) =B n [ A 2n coss n L+B 2n sinB n L+C 2 n cosh Si L+D 2 n sinhB n L 1 


3. System Parameters 


(1) Inertial Moment 

I s =905,443 
1^=6,789,100 
I s 2 =7,086, 601 
I g 3 =145 , 393 
I R *=4,969 
I r 1= 4, 969 
I 2 =9,938 

r 3 


slug- ft 2 
slug-ft 2 
slug-ft 2 
slug-ft 2 
slug- ft 2 
slug-ft 2 
slug-ft 2 
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(2) First Four Modal Coefficients 


Mode No . (n) 

l 

2 

3 

4 

°n 

1.19 

1.29 

1.97 

2.54 

B 'n 

0.033 

0.039 

0.092 

0.152 

w n 

0.274 

0.322 

0.748 

1.24 

A m 

0.161 

0.072 

0.022 

0.068 

B ln 

-0.196 

-0.084 

-0.059 

-0.063 

c m 

-0.168 

-0.075 

-0.023 

-0.068 

D m 

0.195 

0.084 

0.059 

0.063 

A 2n 

-0.039 

0.125 

0.025 

-0.105 

B 2n 

0.069 

-0.196 

0.003 

0.094 

C 2n 

0.058 

-0.167 

-0.025 

0.107 

D 2n 

-0.069 

0.196 

-0.003 

-0.093 

A 3n 

-0.032 

0.003 

0.072 

0.011 

1 II 

1 II 

1 II 

1 II 

! c !! 
1 n> II 

1 CO II 

1 II 

1 II 

I II 

1 II 

1 II 

0. 158E-4 

-0. 109E-5 

-0. 131E-4 

-0. 123E-5 


(3) Other Values 


<Og=7 . 27E- 5 
M=12 . 42 
M R =12 . 42 

X=18. 75 
Y=-32 . 5 
L=- 130 


rad/sec 

slug 

slug 


ft 

ft 

ft 
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, m+'rqrtion - Preliminary Study of th> 

Control Structure^ Ha3S on the Design 

of Control Laws 


_ _ _ + r> 1 1 r + 1 1 r e s are described 

s of large space structui -3 


u s i n g 


where 


The dynamic 
the finite element method as 1 - 


.« • 

M X + CX 

+ KX = BU 

X = nx 1 

vector representing degrees 

M = nxn 

mass matrix 

C - nxn 

damp i ng mat r i x 

K = nxn 

st i f f ness matr i x 

B - nxm 

control influence matrix 

U = m x 1 

control vector 


( 1 ) 


using modal analysis* and modern control theory*, state 
variable feedback control laws of the form 


IJ = - Fp X - FpX 


( 2 ) 


w h e r ■ 


■ r CMS.. , p ar e rate and position control gain matrix 


Fr and F ^ 

of appropriate dimensions are designed. To implement the 
control law given by equation ( 2 ) physical actuator. 
needed. These physical actuators have finite mass and, 
thus, change the mass and stiffness of the structure to b 
controlled. This mass can be as much as fift-^n p^r c«nt 

4 Thus the control laws desi •: 

the uncontrolled structure. 

Without taking this mass into consideration have to be 
reevaluated for their stability and performance degradaf 
AssumingAM and AX are the changes in the mass and stiffn- 


of 

3 ned 
ion. 
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matrices due to actuators the dynamics of the controlled 
system can be written as: 

(M+AM) X + ( C *• B F r ) X + (K + A«+6F p ) X = 0 (3) 

Since the control law is designed for the stability of the 
controlled system, the matrices M, C+BF r , and K are positive 
definite matrices. 5 If the changes in the mass matrix and 
stiffness matrix, AM and AK , are also assumed to be positive 
definite then the matrices (M+AM), ( C + B F r ) , and (K + AK + BF p ) 

are also positive definite. Thus, equation (3) is stable, 
though performance degradation ban not be commented on. The 
assumption that AM and AK are positive is a valid 
assumption, as the dynamics of the oscillatory motion of the 
structure with L he added actuator masses can be described 
based on the finite element method or energy conservation 
techniques, and thus, (M+AM) and (K+AK) must be positive 
definite. As (M +AM), (K+AK) are positive definite and (M+' 
AM), ( K + BFp ) are positive definite, the matrices (M+AM), fc+ 
AK+BFp) are also positive definite. Thus, as M and K are 
positive definite and (M+AM) and (K+AK) are positive 
definite, Am and A K are positive definite. The effect of 
the actuator mass on the structural damping io not 
considered here. 

in this analysis the modal truncation is not taken 
into account and, thus, the control spill -c*r problem will 
not arise. The performance degradation is analysed using a 
two mass, two spring, two actuator system. 
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Numerical Example: 


The two-mass two-spring system is shown in Figure 1 and 
its equations of motion are written as: 



The c o n t r o 1 law of the f o r m 



is designed with the following numerical values for th 
mass, stiffness and control gain matrices. 


m i = 2 , m 2 = !il<l = ^' k 2- 1 

f n = 1, f 1 2 = f 2 1 s 0 ’ f 2 2 = 1 

The numerical simulation is conducted varying the masses and 
stiffnesses, one at a time, and the closed-loop eigenvalues 
are tabulated in Table 1. From Table 1, it can be observed 
that tne change in mass, mo. has a maximum effect on the 
degradation of the closed -loop eigenvalues. A 15% change in 


m 2 
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actuator masses affect the performance of the control law 
that is designed without taking these masses into 
consideration. It is also worthwhile to observe that a 
change in the stiffness moves the eigenvalues to the right 
as well as to the left and can be explained as the effect of 
the increase in the stiffness on one mass or the other. 

To understand the performance degradation due to 
actuators an exhaustive simulation of the closed-loop 
controlled system has to be done with the following 


cons i derat i on 3 : 


1. The masses needed to implement specific control 
rorces have to be evaluate d . 

2. The change in stiffness due to change i n mass 
h as to d e d e t e r m i n e d . 


Simulation h a s to be condu c t e d w i t h 
total mass and st if mess matri ces r 


hanges in the 


h e r t h a n 


i n d i v i d u a 
a co n t r o 1 s y s 
the actuator' masse 
i ncorporat i ng the 
d y n a m i c model u n t i 


1 ma sees as is done i n t h \ s s t u J y • 
t e m ci t* c i g n to acco rn m o d a t e t h e effect of 
s has to be done in iterati ve fashion 
c h a n g e i n m a s s a n d s t i f f n e : ; s i n t o t h e 
I a s a t i s f a c f o r' y o o n t r o ! i -a w i 3 a r' r i v e d 


at 


1 . 


Z i e n '• i e v / ■ 
i ] I do C'd 
M .a i" ! , •• r < 1 


R e f e r e n o a s 

r < o.C.. I ' i ^ fi nite E i erne nt Method. • i'icGraw- 

C .. 'Hip on v Inc . . N — w V err . I ■> > 

• J . . Genera.! .! JLLJ. i±S±Jliid..a *!. jOiLLILLiJin AO , rf , 1 i aL , 
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K wake r n a a k < J ■ and s i v 3 n ^ R . • .L»J — n..e...™.j — im.. 
S ystems , John Wiiey & Sons. New York. 1370. 
£ 3 r' t h Pointing Satellite ( E Pi ) s 1 1 ' u c t u 1 « 0 1 
NASA Internal Document, Jan.. 1989. 

Bellman, R . , IiijLLsd^ M a t rix AnaLa Jss 

Hill Book Company Inc., New York, 1960. 


Q-l. 


scr i pt i on , 


i 3 . M c G i" a w 
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Figure 1: Two Mass-Two Spring System 
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Percentage change in 
masses and stiffnesses 


Closed loop eigenvalues 
(complex conjugate pairs) 


ra l 

m 2 

k l 

k 2 

1 



2 



0 

0 

0 

0 

-0.472 

+ 

jO . 710 

-.277 

+ 

jO . 163 

15 

0 

0 

0 

-0 . 464 

+ 

jO. 709 

-.254 

+ 

j0.154 

0 

15 

0 

0 

-0.420 

+ 

jO. 684 

-.268 

+ 

jO . 163 

0 

0 

15 

0 

-0.477 

+ 

jO. 729 

-.272 

+ 

jO. 171 

0 

0 

0 

15 

-0.465 

+ 

jO. 756 

- .284 

+ 

jO. 168 


Table 1 : Closed-loop Eigenvalues due to Changes in Mass 

and Stiffness Values. 
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CONCLUSIONS AND RECOMMENDATIONS 


V . 

The maximum principle of Pontryagin has been applied to 
the rapid maneuvering problem of the planar, flexible 
orbiting SCOLE. The resulting two-point boundary value 
problem is solved by applying the quasi 1 i neari zati on 
technique, and the near-minimum time is obtained by 
shortening the maneuvering time in a sequential manner until 
the controls are near the bang-bang type. The results 
indicate that responses of the nonlinear system for the 
flexible modal amplitudes may be significantly different 
from those of the correspond!' ng linearized system for rapid 
slewing maneuvers. This research is currently being 
extended to the three dimensional slewing of the flexible 
SCOLE system. 

From an analysis and simulation of the SCOLE station- 
keeping dynamics it is found that the flexible vibrations of 
the mast are not greatly excited during typical station- 
keeping operations. System responses are highly dependent 
on the force actuator locations and the numerical values of 
the state and control penalty matrices included in the LQR 
control law design. Force actuators mounted at 1/3 and ~/3 
of the mast length along the mast are effective in 
suppressing the flexible mast viorations. 

A preliminary examination of the effect of actuator 
mass on the design of control laws for large flexible space 
3 y s t e rn s d e mo n s t r a t e s t h a t act u a t o i m asses can influence the 
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performance of the closed-loop system where the control law 
has been designed without taking these masses into 
consideration. To understand better the possible degrada- 
tion in performance due to the presence of actuator masses 
additional studies are required to accurately evaluate the 
changes in the stiffness matrix due to specific actuator- 
masses, and simulations must be performed incorporating 
changes in the total mass and stiffness matrices, rather 
than individual masses as was done here. 

Finally, the current (1989-90) grant work has been 
redirected so as to lend greater support to the new 
Control s/Structures Interaction (CSI) program and focusing 
on specific CSI evolutionary conf i gurat i ons , in addition to 
the treatment of the SCOLE 3-D slewing problem. 
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